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Abstract 

In recent years there has been an explosive development of interest in 
the measurement of forces at the microscopic level, such as within living 
cells 013 El, as well as the properties of fluids and suspensions on this 
scale 4 , using optically trapped particles as probes. The next step would 
be to measure torques and associated rotational motion 0- This would 
allow measurement on very small scales since no translational motion 
is needed. It could also provide an absolute measurement of the forces 
holding a stationary non-rotating particle in place. The laser-induced 
torque acting on an optically trapped microscopic birefringent particle 
[Sj can be used for these measurements. Here we present a new method 
for simple, robust, accurate, simultaneous measurement of the rotation 
speed of a laser trapped birefringent particle, and the optical torque acting 
on it, by measuring the change in angular momentum of the light from 
passing through the particle. This method does not depend on the size 
or shape of the particle or the laser beam geometry, nor does it depend 
on the properties of the surrounding medium. This could allow accurate 
measurement of viscosity on a microscopic scale. 

1 Introduction 

Optical torques have been measured previously; two methods have been used. 
Firstly, if a particle with known birefringent properties has a simple and 
accurately known size and shape, ideally a flat disc, the torque can be calculated 
from the beam power jjjj. However, particles of more complex shapes will 
often be used in experiments or encountered in samples; for example, spherical 
particles are ideal for making measurements of viscosity. Secondly, torques 
have been determined by measuring the rotation speed in a medium of known 
viscosity @G]. This method cannot be used if the aim is to measure an unknown 
viscosity. This method would also fail if there were other torques acting on the 
particle, or if the viscous drag is affected by nearby walls or other particles, 
or if the particle is not rotating. Previous methods for measuring rotation 
speeds, based on the periodic variation of backscattered light |8], can also have 
problems. Very regular or rotationally symmetric particles provide insufficient 
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variation or variation at an increased frequency. Consideration of the basic 
physical processes giving rise to the torque gives a new method for measuring 
the torque and rotation speed that overcomes all of these problems. 

Since optical torques and forces are very small, microscopic particles are 
ideal for the observation and application of optical torques and rotation. Such 
microscopic particles will typically be confined within a laser trap. Strongly 
focussed laser light incident on a transparent particle, usually in a liquid 
medium, will produce a gradient force acting on the particle towards the region 
of highest irradiance. If this gradient force near the focus is stronger than 
scattering and absorption forces, the particle will be trapped at the beam 
focus, where the irradiance is highest. This technique of three-dimensional 
confinement and manipulation is called laser micro-manipulation, trapping, or 
optical tweezers. 

2 Polarised beams 

A monochromatic laser beam can be written as a plane wave in terms of two 
orthogonal components: 

E = (E x x + E y y) cxp(ikz — iwt) (1) 

where the beam is propagating in the z-direction. In general, the amplitudes 
E x and E y are complex in order to account for the phases of the components. 

There are two cases of special interest. The first is when the phase angle 
between the complex amplitudes E x and E y is equal to or it , in which case 
the total electric field always lies in a single plane resulting in a beam which is 
linearly polarised. The direction of the x-axis can be chosen to coincide with 
the plane of polarisation, so the beam can be written as 

E = E p xexp(ikz — iu>t) (2) 

where E p is the complex amplitude of the linearly polarised light 

The second special case is when the phase angle is ±7r/2, and \E X \ = \E y \. 
In this case, E y — E x expi9, with 8 = ±tt/2. The total electric field has a 
constant magnitude, with the direction varying with the optical frequency uj so 
that the beam is circularly polarised. When 9 = +tt/2, the electric field has 
a positive, or right-handed, helicity. Such a beam is here called left circularly 
polarised. When 9 = —it/2, the beam has negative helicity and is called right 
circularly polarised. A circularly polarised beam can always be written as 

E = (E c x ± iE c y) exp(ikz — iu)t) (3) 

with the sign depending on whether the beam is left or right circularly polarised, 
and E c is the complex amplitude. 

In general, however, the phase angle 9 will have a value between these 
limiting values, or even if 9 = ±7r/2, \E X \ ^ \E y \. In these cases, the beam is 
clliptically polarised, and the electric field vector E traces out an ellipse during 
each optical period. 

Recognising that we can rewrite equation © for a circularly polarised beam 

as 

E = E c (x ± iy) exp(ikz — iut), (4) 
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El = —(E x -iE y ) 



we see that any beam can be represented as a sum of two circularly polarised 
components using the (complex) orthogonal basis vectors 

e L = -^=(x + «y) 

e« = -j=(x-iy) (5) 

as 

E = (E L e L + E R e R ) exp(ikz — iwt). (6) 

The amplitudes of the left and right circular components can be found from the 
x and y amplitudes in the linear orthogonal representation (equation Q): 

1 

Er = ^(E x + iE v ) (7) 

When \El\ = \E R \, the beam is linearly polarised, with the plane of polari- 
sation given by the phase angle between the complex amplitudes El and E R . If 
El = (and E R ^ 0), the beam is right circularly polarised, and left circularly 
polarised if Er = 0. 

The time-averaged irradiance is given by 9 

ce E* L E L ce E R E R 
1 ~ ^ 1 2 = 1 L + 1 R- (p) 

Although we have only considered the beam as a classical EM wave so far, the 
fact that the angular momentum of left and right circularly polarised photons 
is ±h can be used to simply find the angular momentum of the beam. Since the 
energy of a photon is Tiuj 1 the photon flux per unit area is 

flUJ flUJ flUl ' ^ 

giving an angular momentum flux per unit area of 

L z = (I L - Ir)/u. (10) 

Thus, the beam can be considered to have a net circularly polarised component 
with a power of \Il — Ir\ which contributes to the angular momentum of the 
beam, and a linearly polarised component of I — \Il — I R \ = 2 min(/i, I R ) which 
does not contribute to the angular momentum of the beam. We can define a 
coefficient of circular polarisation a z by 

o z = (I L - Ir)/I, (11) 

and write the angular momentum flux density of the beam as 

L z = aj/u). (12) 

When the irradiance is integrated across the whole beam, the total power 
can be obtained and will be given by P = Pl + Pr = J IidA + J IrcIA. A 
suitable average coefficient of circular polarisation can be defined by 

a z = (P L - Pr)/P, (13) 

with the resulting total angular momentum flux of the beam being 

L z = a z P/uj. (14) 
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3 Optical torque 

If the beam passes through some birefringent material, the polarisation will be 
affected. In general, <r z will change. The incident beam will have an initial 
coefficient of circular polarisation a Z i n , and will have an emergent polarisation 
described by o~ zout . Thus the angular momentum of the beam will change, and 
a reaction torque on the birefringent material will result. The reaction torque 
is equal to the change in the angular momentum flux: 

T = {o zln - <7 zout )P/u! (15) 

assuming that absorption and reflection can be ignored. Equation (|15|l is 
general, and can always be used to find the torque if the coefficients of circular 
polarisation of the incident and outgoing beams are known or can be found. 

Although equation l|15|l applies in general, it is instructive to carry through 
a detailed calculation for a simple case: a uniform sheet of birefringent material, 
for example calcite. 

A uniaxial birefringent material such as calcite can be described by two 
refractive indices: an ordinary refractive index n Q for electric fields normal to 
the optic axis, and an extraordinary refractive index n e for electric fields parallel 
to the optic axis. For calcite, n — 1.66 and n e — 1.49. Consider a thickness d 
of uniaxial birefringent material with the optic axis in the xy plane, at an angle 
of 9 to the cc-axis. The front face of the material is at z — zq, and the rear face 
is at z = zq + d. If the electric field of the incident beam at the front surface 
of the material is given by equation Q), we can express this in terms of unit 
vectors i and j parallel to and normal to the optic axis: 

E = \(E X cos 9 + E y sin 9)\ + {-E x sin 9 + E y cos 6>)j] 

x exp(i/czo — iujt). (16) 

In terms of circular components, this gives 

E = -^=[{E X - iE y ) exp(i6)e' L + (E x + iE y ) cxp(-i9)e' R ] 
v 2 

x exp(ikza — iut) (17) 

where e' L = ^(i+y) an d e' R = -i=(i — ij) The coefficient of circular polarisation 
is given by 

= E* L E L - E* R E R 

(18) 



i(E x E* y - E*E y ) 



E%E X + E*E y 

After passing through the thickness d, the field will be 

E = [{E x cos 9 + E y sin (9) exp(ikdn e )i 

+ (—E x sin 9 + E y cos 9) exp(ikdn a )j] exp(ikzo — iujt), (19) 

which we can express in terms of circular components 

El = -^=[(E X cos 9 + E y sin 9) exp(ikdn e ) 
v2 
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—i(—E x cos9 + E y sin#) exp(ikdn )] 

Er = —=[{E X cos 9 + E y sin 9) exp(ikdn e ) 
V2 

+i(—E x cos 9 + E y sin 9) exp(i/cdn )] (20) 

We define the convenient notation 

A = kd(n - n e ) (21) 

The coefficient of circular polarisation of the emergent light is 

(T Z out = [i cos A(E X E* — E*E y ) 

- sin A{(E*E X - E*E y ) sin 29 - (E X E* + E*E y ) cos 20}] 
/(E*E x +E* y Ey) (22) 

giving a torque per unit area of 

r = ^[i(E x E;-EZE v )(l-cosA) 

+ sin A{(E*E X - E*Ey) sin 29 - (E X E* + E*E y ) cos 29}} (23) 

If the incident light is linearly polarised (E y = 0), the torque is 

r = — sin AE*E sin 29, (24) 
2uj 

which acts to align the slow axis of the particle with the plane of polarisation if 
n > n e , or normal to the plane of polarisation if n e > n . If the incident light 
is left circularly polarised (E y — iE x ), the torque is 

t = — E%E (l-coaA) (25) 
to 

which is independent of the orientation of the birefringent material. 

If the birefringent material is of a uniform thickness, the total torque can 
be simply calculated from this (Hj- In general, a laser trapped birefringent 
particle will have a varying thickness, and direct calculation of the torque will 
not be feasible. Also, if the orientation of the birefringent particle is different, 
so the the optic axis does not lie in the xy plane, the calculation will be further 
complicated. Equation (|15(l . however, is general, and will still apply, and the 
torque acting on the particle can be deduced from the change in polarisation of 
the light. 



4 Optical torque measurement 

Consider a circularly polarised laser beam used to trap a microscopic particle 
composed of a uniaxial birefringent material such as calcitc or a suitable 
polymer. If the optical torque is large enough to overcome forces holding the 
particle in place, the particle will rotate at a speed determined by the equilibrium 
between the optical torque and other forces such as viscous drag. In this way, 
a probe particle can be used to measure viscosity on a microscopic scale. If the 
particle does not rotate, the optical torque can be used to determine the torque 
due to static forces acting on the particle. 
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The maximum torque and rotation rate will occur when the incident beam 
is completely circularly polarised (ie u Z i n = 1)- The torque in this case will also 
be constant as well as maximal 6 , and we will only consider this case here. 
The torque r acting on the trapped particle is given in equation (|15|) by the 
difference between the incident and outgoing angular momentum fluxes, and in 
this case, assuming no reflection or absorption, is 

t = (1 - a zout )P/uj. (26) 

Measurement of the outgoing polarisation o zout and beam power P gives an 
absolute measurement of the torque, which does not depend on the mechanical 
properties of the surrounding medium or the particular size or shape of the 
particle or laser beam. 

We can also note that the plane of polarisation of the linearly polarised 
component of the outgoing beam exiting a rotating birefringent particle will be 
rotating at the same rotation rate O as the particle. If the outgoing beam is a 
(rotating) purely plane-polarised beam, as would occur if the particle acted as 
a quarter-wave plate, rotating at O, and of power P, and is passed through a 
linear polariser, the measured power P m will be P m = (1 + cos20i)P/2 (with 
variation at a frequency of 20 since a rotation of 180° rotates the plane of 
polarisation onto itself) ^H]- By measuring this power, the rotation rate O of 
the trapped particle can be simply determined. This will still be the case for an 
elliptically polarised beam, as the same variation at a frequency of 20 will be 
observed. The angular momentum associated with this rotation of the plane of 
polarisation will be negligible asl]<Cw. 

In the general case, there will be an elliptically polarised outgoing beam, 
consisting of both plane and circularly polarised components. The measured 
power P m after the outgoing beam passes through a linear polariser acting as 
an analyser will be 

P m - {1 + (1 - * 2 zout ) l/2 cos20i}P/2. (27) 

Measurement of the variation of the transmitted power therefore allows the 
determination of the rotation period of the trapped particle, and the degree of 
(but not the direction of) circular polarisation. The result of a measurement of 
this type will be as shown in figure 1. This measurement is an average over the 
beam, and it is not important whether or not the entire beam passes through the 
particle. In the case where the particle is not rotating, due to some restraining 
torque, the plane of polarisation of the transmitted light will not be rotating. 
The degree of circular polarisation can be measured in this case by rotating the 
linear polariser, which will give the same result where O is the rotation rate of 
the polariser relative to the particle. The orientation of the particle can also be 
determined from the position of the measured power maxima. 

In many cases, the direction of the transmitted polarisation will be known 
beforehand — such as when the particle is insufficiently thick to change the 
direction of polarisation (note that a calcite particle approximately 3fim thick 
is a A/2 plate for 1064 nm light), or when the particle is small and the outgoing 
light is dominated by light that did not pass through the particle and has not 
changed in polarisation. If necessary, the direction of circular polarisation can 
be measured simply by placing a reversed circular polariser (eg a quarter-wave 
plate followed by a linear polariser appropriately oriented) in the beam path 
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Figure 1: The power which would be measured through a plane polariser after 
the beam has passed through a birefringent particle is shown for a beam with 
&zin = +1, &zout = +0.7 and particle rotation frequency SI — 10 Hz. The 
mean power measured through the polariser is half of the power incident on the 
particle. The frequency of the variation is two times the rotation rate of the 
particle. The optical torque can be found from the amplitude of the variation 
and the measured power once the direction of the transmitted polarisation is 
known. In this case, for a 100 mW trapping beam of wavelength 1,064 nm, an 
optical torque of 16.9 pN-/im is being exerted. 



instead of a linear polariser. In the case where the trapping beam has a right- 
handed helicity (left circularly polarised, with cr Z i n = +1), the light emergent 
from the particle can be described in terms of left and right circularly polarised 
components Pl and Pr, where Pl = (1 + u zou t)Pj2 and Pr = (1 — a zout )P/2. 
If the output beam is predominantly left circularly polarised, a ZO ut > 0, and 
Pl > Pr- A right circularly polarised beam has a zou t < 0, and Pr > Pl- It is 
only necessary to determine which of these two components is larger, rather than 
to measure each one individually since \<J zou t\ is already known. In this way, 
the direction of circular polarisation can be determined, and a zou t as opposed 
to merely \a zou t\ can be found. Once a zout is known, the optical torque acting 
on the particle can be found using equation 1)150. A measurement of this type 
will be as shown in figure 2. 

It should be noted that this technique is robust. It is not necessary 
to measure the power of the entire transmitted beam; it is sufficient to 
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Figure 2: The two circularly polarised components of the transmitteed light 
can be measured to determine the direction of circular polarisation. Pl is the 
power of the left circularly polarised component, and Pr is the power of the 
right circularly polarised component. 

measure the portion of the beam that has passed through the trapped particle. 
Similarly, reflections are not likely to cause significant error. Some of the 
incident beam will be reflected from the trapped particle; the reflection 
will depend on the angle of incidence and the refractive indices of the 
particle and the surrounding medium. For example, for calcite trapped in 
water, the Fresnel amplitude coefficients for reflection at normal incidence are 
(n wa ter - n ca i c ite) I (n wa ter + n ca i cite ), which gives reflected amplitudes of 
—0.057-2^ and —0.11iE x for linearly polarised components parallel to and 
normal to the optic axis respectively. In terms of circular components, this 
becomes El = — O.Q2Elq and Er — —0.08Elo, showing that the torque due to 
backreflected light will be less than 0.6% of the available torque. Therefore, the 
reflected light will not cause any significant error. 

5 Conclusion 

A simple method of measuring the rotation speed and the optical torque applied 
to a laser trapped birefringent particle has been described. This method can 
be used even if the viscosity of the medium in which trapping is performed 
is unknown, and provides a means to measure this viscosity. Thus, this 
method is suitable for employment in a micro-rheometer, which could be simply 
constructed by trapping a birefringent probe particle in the fluid of interest. A 
suitable test particle would be a small fragment of calcite, the exact shape not 
being critical at the very low Reynolds numbers encountered in these cases, or 
a more ideal shape could be fabricated from a birefringent polymer JJ. As 
the optical torque can be controlled by varying the power, the probe particle 
rotation speed can be varied, allowing, for example, the investigation of non- 
linear properties of the fluid. Wall effects and the small-scale behaviour of 
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polymer and colloidal suspensions could be investigated, or even the rheological 
properties of intracellular fluids or membranes in vivo. 
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